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1 Introduction 



Gaudin type models [T] have important applications in physics. For instance, the XXZ 
Gaudin model has played an essential role in the study of the reduced BCS model whose 
exact solutions were first found by Richardson [2]. In fact, the conserved operators for the 
BCS model can be mapped to a set of XXX Gaudin Hamiltonians in a non-uniform magnetic 
field and the BCS Hamiltonian is expressible in terms of these operators [HJHUH]. This result 
also exposes a relationship between the BCS Hamiltonian and the perturbed WZNW model 
at critical level [6] based on the connection of Gaudin models with the solutions of the KZ 
equations [7J E\- 

Eigenstates and the corresponding eigenvalues for the open XXZ Gaudin model for the 
boundary conditions with three free boundary parameters were derived in Ref.jH]. In this 
paper, we consider the most generic boundary conditions specified by the non-diagonal K- 
matrices in [TUl [EE] , leading to the open XXZ Gaudin model in this paper which depends on 
four free boundary parameters. We compute the scalar products of this Gaudin model, and 
give their explicit expressions in terms of determinants, that is the determinant representa- 
tions of the scalar products. 

This paper is organized as follows. Section 2 provides some preliminaries on the boundary 
inverse scattering method. In section 3, we briefly describe the open XXZ Gaudin magnet 
associated with non-diagonal boundary K-matrices. In section 4, we derive the Bethe ansatz 
equations for the open XXZ Gaudin model and the symmetric, polarization-free expressions 
for the pseudo-particle creation operators. In section 5, we obtain the determinant represen- 
tations for the partition and correlation functions of the model. We summarize our results 
in section 6 and present the details of some derivations and proofs in the Appendices. 

2 Preliminaries: the inhomogeneous spin-| XXZ open 
chain 

Let V be a two-dimensional linear space and cr ± , a z be the Pauli matrices which give the 
spin- 1 representation of su{2) on V. The spin-^ XXZ chain can be constructed from the 
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well-known six-vertex model R-matrix R(u) G End(V <8> V) [12] given by 

(2.i; 



R(u) 



( 1 \ 

sintt sinr? 

sin(u+r]) sin(«+?j) 

sinr? sintt 

sin(u+r?) s'm(u+r)) 

V 1 / 



Here we assume rj is a generic complex number. The R-matrix satisfies the quantum Yang- 
Baxter equation (QYBE), 

Ri j2 (ui - u 2 )Ri,3(u 1 - u 3 )R 2i3 (u 2 - u 3 ) = R 2 , 3 (u 2 - u 3 )R li3 (u 1 - u 3 )R lj2 (u 1 - u 2 ), (2.2) 

and the unitarity, crossing-unitarity and quasi-classical properties [9]. We adopt the standard 
notations: for any matrix A G End(V) , Aj is an embedding operator in the tensor space 
V <g) V <g) • • -, which acts as A on the j-th space and as identity on the other factor spaces; 
Rij(u) is an embedding operator of R-matrix in the tensor space, which acts as identity on 
the factor spaces except for the z-th and j-th ones. 

One introduces the "row-to-row" (or one-row) monodromy matrix T(u), which is an 2 x 2 
matrix with elements being operators acting on V® N , where iV = 2M (M being a positive 
integer) , 

T (u) = R , N (u - z^Rq^^u - zat-i) • • • Rq,\{u - z x ). (2.3) 

Here {zj\j = 1, • • • , N} are arbitrary free complex parameters which are usually called inho- 
mogeneous parameters. 

Integrable open chain can be constructed as follows [32]. Let us introduce a pair of 
K- matrices K~(u) and K + (u). The former satisfies the reflection equation (RE) 

#1,2(1*1 - u 2 )K{(ui)R 2)1 (ui + u 2 )K 2 (u 2 ) 

= K^{u 2 )R 1)2 { Ul + u 2 )Ki(u 1 )R 2 , 1 (u 1 - u 2 ), (2.4) 

and the latter satisfies the dual RE 

Ri, 2 (u 2 - ux)K^ (u 1 )R 2il (-u 1 -u 2 - 2rj)K%{u 2 ) 

= K+{u 2 )R 1>2 {- Ul -u 2 - 2r t )Kt{u l )R 2 , l {u 2 - Ul ). (2.5) 

For open spin-chains, instead of the standard "row-to-row" monodromy matrix T(u) ( 12. 3p . 
one needs to consider the "double-row" monodromy matrix T(u) 

T(u) = T(u)K~(u)f(u), f(u) = r _1 (-«). (2.6) 



Then the double-row transfer matrix of the XXZ chain with open boundary (or the open 
XXZ chain) is given by 

t(u) = tr(K + (u)T(u)). (2.7) 

The QYBE and (dual) REs lead to that the transfer matrices with different spectral param- 
eters commute with each other |13j : [t(u),t(v)\ = 0. This ensures the integrability of the 
open XXZ chain. 

3 XXZ Gaudin model with generic boundaries 

We will consider the K-matrix K~(u) which is a generic solution to the RE ( 12 .4p associated 
the six-vertex model R-matrix [101 E] 

K-(u)=( k M %"l)=K{u). (3.1) 



k\{u) //>; a 

The matrix elements are 



k\{u) 
k\{u) 
kf(u) 



cos(Ai — A 2 ) — cos(Ai + A 2 + 2£)e" 
2 sin(Ai + f + u) sin(A 2 + f + u) 
-ism(2u)e~ i(Xl+X2) e- iu 
2 sin(Ai + f + u) sin(A 2 + f + u) ' 

ism(2u)e iiXl+X2) e~ iu 
2 sin(Ai + I + u) sin(A 2 + \ + u) ' 



-2iu 



k 2 (u) = cos(A! - A 2 )e- 2 ^ - cos(A 1 + A 2 + 2fl 
2[U) 2sin(Ai + e + M)sin(A 2 + e + M) ' 
The corresponding dual K-matrix K + (u) is a generic solution to the dual reflection equation 
(12. 5p with a particular choice of the free boundary parameters: 

k + \(u) k+l(u) 
k + l(u) k+ 2 2 (u) 



with matrix elements 



k+ i,. _ cos(Ai - \ 2 )e- ir > - cos(Ai + A 2 + 2Qe 2 ^ 
2 sin(Ai + £ — w — 77) sin(A 2 + £ — it — 77) 
i sin(2w + 2r/)e- i(Al+A2) e i "- i ' 7 



2, 



fc + i(iz 



2 sin(Ai + £ — w — 77) sin(A 2 + £ — m — 77) 

-isin(2ti + 2^?)e ^(Al+A2) e ^ " +^, ' 
2 sin(Ax + £ — u — 77) sin(A 2 + £ — u — 77) 



. +2 , , cos(A 1 - \ 2 )e 2iu+ir > - cos(A! + A 2 + 2£) e - i " 
2 sm(Ai + £ — it — 77) sm(A 2 + £ — w — 77) 
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The K-matrices depend on four free boundary parameters {Ai, A2, £, £} which specify inte- 
grable boundary conditions [11]. We remark that K~(u) does not depend on the crossing 
parameter r] but K + (u) does. The parameter £ is required to have the following expansion: 

S = £ + T,A + 0(ri 2 ), rj^O. (3.5) 

This results in the relation, 

\im{K + (u)K-(u)} = \im{K + (u)}K(u) = id. (3.6) 

Let us introduce the generalized XXZ Gaudin operators [lj {Hj\j = 1, 2, ■ • • , N} associated 
with the spin-| XXZ model with boundaries specified by the K-matrices fl3.ll) and (13. 3ft : 

2M 



1 f a z a z -1) 

Hi = Tfa) + E _ ^ + a k a+ + cos( Zj - z k )-^— j 

where ^-(u) = ^{/^(m)}^^), j = l,---,iV, with = tr {K + (M)i? 0j (2n)F 0j }, 

and {zj} correspond to the inhomogeneous parameters of the spin- ^ XXZ chain with generic 
open boundaries. For a generic choice of the boundary parameters {Ai, A 2 , £, £}, r 3 (w) is an 
non-diagonal matrix, in contrast to the situation in [3J. 

The XXZ Gaudin operators ( 13. 7p can be obtained by expanding the double-row transfer 
matrix r{u) (12. 7p at u — Zj around r] = [9]: 

T( % ) = id + ^ + 0(r ? 2 ), j = l,---,N, (3.8) 
^ = ^ r ^)U- ( 3 - 9 ) 
Then the commutativity of the transfer matrices {t(^)} for a generic r\ implies that 

[H j} H k ] = 0, i,j = l,---,N. (3.10) 

Thus the Gaudin model with the local Hamiltonian ( 13. 7\\ is integrable. 

4 Eigenstates and the corresponding eigenvalues 

The relation ( 13. 9 j) between {-Hj} and {r(z,,)} and the fact that the first term of (13. 8p is 
the identity operator enable us to extract the eigenstates of the Gaudin operators and the 
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corresponding eigenvalues from those of the XXZ chain obtained in [HI [151 El HI CEO EES [HI 
US [211 E3 [23]. 

Let us introduce the states |f^) and |^ 2 )), 

l^ (1) ) = ( ! J®---®^ 6 j J, (4.1) 
l^ (2) )= j ®---®( e x , (4.2) 



and the matrix g(-u) G End(F) and the associated gauged Pauli operator a ± {u) G End(V) 

/ -i(u+2A0 -i(«+2A 2 ) \ 
<K«) = ( 1 I J' ( 4 - 3 ) 

^(u) = giu^giuY 1 . (4.4) 

Then we define the states, 

M 

i{« i (i) }> w =n s ( v * (i) )i n(i> >. ( 4 - s ) 



M 



iK (2) }) (2) =n^ (2) )i fi(2) )- ( 4 - 6 ) 



i=l 



The associated operators B{u) and C{u) are 

= V , ^(A 1+ e + ^)sin(A 2 + e- ,^sin(2,) 

^— j' sm(Ai + £ — wj sm(A 2 + q — u) sm(u — sin(« + £j) 

^— j' sm(Ai + £ + u) sm(A 2 + c, + u) sin(u — sm(« + Zj) 
Using the same method as in [24J, we can show that the above states (14. 5 p and (14. 6p are 
the common eigenstates of the Gaudin operators {Hj} given by (13. 7p provided that the 

parameters {i^} satisfy the following two sets of Bethe ansatz equations 

1 - A 1 + A 

sin(Ai + f + v { a ] ) sin(Ai + £ - wi 1} ) sin(A 2 + f + wi^) sin(A 2 + £ - u^) 

AT 2 2M ^ 

= • / (i) W\ ■ i (i) , (T)T — 5^ • / (!) \ . / (i) i r> ( 4 - 9 ) 

fc ^ Q sm(t£ - ffc ) sm(fA, + ^ J fc=1 sm(t>^ - z fe ) sin(% + z&) 

1 + A 1 - A 

sin(A! + i + vi 2) ) sin(Ai + £ - v^) sin(A 2 + f + wi 2) ) sin(A 2 + f - vi 2) ) 

M 2 2M ^ 



• / (2) (2)n . / (2) . (2), • / ' / 

k ^ a sm(t£ - ffc ) sm(t>^ + v k J ) k=1 sm{Va - z k ) sm(% + z&J 

a = l,---,M. (4.10) 
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Here A is the parameter of first order expansion of £ in terms of rj, as defined in (13. 5p . 
Namely, 



^ |{«J?»« = Ef |{««}>«, i = l,2, (4.11) 



where are given by 



2 

E Asin(2z,) 
C ° t(Aj + £ ~ %) " s,n(A 1+ C-3 j )sm(A 1+? + % ) 

M • /<") \ 

k=i sm(^ - sm(^ + z,-) 



2 

^C2~i * x Asin(2z,- 
£ 2j = COt2^+> cot (Aj- + Z-Zj)-— -r rVj 



At 



sin(A 2 + £ - zj) sin(A 2 + £ + Zj' 



Esin(2^ 
m , (2) • ( 4 - 13 ) 

k =i sm K - z i) sm K + h) 
5 Determinant representations of the scalar products 

To obtain correlation functions, it suffices to calculate the scalar products of on-shell Bethe 
states with general off-shell Bethe states [12] (see also [251 EE] for the open XXZ chain with 
diagonal boundaries). In this section, we will obtain the explicit expressions of the following 
scalar products for the open XXZ Gaudin model with non-diagonal boundary terms: 

^• 2 (K}; K (2) }) = (1) <K}I-N 2) }> (2) , s 2,1 (K}; = (2) <K}IK W }) (1) , (5.i) 

^(K); K (1) }) = «(K}|K W }}«, S^({u Q }; {v? }) = V{{u a }\{vf } }) (2) , (5-2) 
where ( - 1 - ) ({'Uq}| and ^ are defined by 

(1) (K}| = <n (1) |C(ti Ar )...G(tt 1 ) l (5.3) 

(2) (K}| = (Q {2) \B(u M ) . ..B{ Ul ), (5.4) 

with being the dual states of |fi (2) ), respectively, 

(tf>|= (1, -e-^ 1+2A2 ) ) ( 1. -e-<(**+2*0 ) , (5.5) 

' 2sin(Ai - A 2 ) I 



^ • -i( 2j +Ai+A 2 ) 



(^( 2 )| = m^_ — W-l, e-^+ 2Al ) )«)••• ® ( -1, e-^ +2Al ) ) . (5.6) 

I - L j L 2 sm(Ai — A 2 ' 1 



By means of (I4.1l) - (j4.8p . ( 15. 5 p and (15. 6p . the scalar products can be written as 

S 1,a ({ua}; M) = (f\C(u 1 )---C(u M )C(v 1 )---C{v M )\^), (5.7) 

^(WlM) = {^\B{u 1 )---B{u m )B{v 1 )---B{vm)\ t), (5.8) 

^(K};^}) = (t|C , (Mi)-"^(«M)S(t;i)-- -5(^)1 1r>, (5-9) 

S 2,2 (K};M) = (^| j b(«i)---S(wm)^i)---^m)|^>. (5.10) 

Here | ff) and (ff | (resp. | JJA and (Jj |) are the all spin-up state and its dual (resp. all 
spin-down and its dual), and C(u) and B{u) are given by 

~ ^ sin(Ai + £ + sin(A 2 + £ - Zi) sin(2w) _ 

#(«) = > — 7T t , . , , . . , . . , r X (J- , (5.11) 

' sm(Ai -\- 1, — u) sm(A 2 + q — u) sm{u — Zi) sm{u + Zi) 

C{u) = V , " ^(Ai + $-^)sin(A 2 + e + ^)sin(2 M ) ^ x 

^ sm(Ai + £ + u) sm(A 2 + £ + u ) sm(u - sm(u + Zi) 

5.1 Scalar products S 1,2 and S* 21 

Let us introduce two functions 

Z$({uj}) = S 1 ' 2 ({u a };{v l }) = (H\C(u 1 )---C(u M )C(v 1 )---C(v M )\^ 1 (5.13) 
Z§\{uj}) = S 2 '\{u a };{v i }) = (ty\B(u l )---B(u M )B(v l )---B(v M )\1t), (5.14) 

where the N parameters {u~j\ J = 1, . . . iV} are defined as 

Ui — Ui for i — 1, . . . M, and %+i = «i for z = 1, . . . M. (5.15) 

Note that these functions Z^\{uj}) and Z^\{uj}) may correspond to the partition func- 
tions of the Gaudin model with domain wall boundary condition and one reflecting end 
[27J EE] specified by the non-diagonal K-matrices (13.11) and (13. 3p . 

We find that the functions Z^ \{uj}) above can be expressed in term of the determinants 
of the N x N matrices jVW({u a }] {zi}) a j, 

ZP({u }) = ^=1^1^^ + Zj)sm(u a - Zi)det^ k \{u a }]{zi}) a j ^ ^ 



n a>/3 sin(u a - up) sin(w Q + up) Y[ k<j sin(z k - zj) sin(z k + 



1) 



where 



yV>(Ui\-(rX\ sintAx +Z-Zj) sin(A 2 + £ + z t ) sin(2^) 

<sr \ z ij)a,j ~ . / > I _x . 77 , t , -\ • 2/- \ • 2/- , \' l - 1 '/ 1 

sm(Ai + 4 + «) sm(A 2 + £ + u) sin (it — Zi) sin (it + ZjJ 

^»({«.}; WW = = 7T ;"'' A i^^-)7( A ;^-^) 3i f S ) (5.18) 

sm(Ai + £ — u) sm(A 2 + £ — u) sin [u — Zi) sin (u + Zi) 

In appendix A we give the proof of this determinant representation of the partition functions. 
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5.2 Scalar products S 1,1 and S 2,2 

In this subsection, we calculate the scalar products S 1 ' 1 ({u a }; {v^}) (15. 9p and S 2,2 ({u a }; {v^}) 
(I5.10p . In this case, we assume that {v^} satisfy the associated Bethe ansatz equations. 

Let us first consider the scalar product {/uj }) ( 15. 9p . We insert in the scalar 

product a sum over the complete set of states \ji, • • ■ , j m ) between each operator, where 

' • " is the state with i spins down at the sites ji, • • ■ ,ji and 2M — i spins up at the 
other sites. We are thus led to considering the intermediate functions 

G w Oi, • • • , Ui\j i+ i, j M ; = (ji+i, ■ ■ -,jM\C{ui) ■ ■ ■ C{ux) 

xB(vi 1) )...B(v^ ) )\f), (5.19) 

which satisfy the following recursive relations: 

G^(u u ---,u l \j l+u ---,j M ;{vl 1) }) 
= • " " ihi\C(ui)\j,j i+u - ■ ■ ,j M ) 

x G^im, • • • , Ui-xljJi+u ■ ■ -,jm; {v^}), z = 1, • • • , M. (5.20) 
The last one of these functions is the scalar product ^({ita}; {v\ }), namely, 

G (M) («i, • ■ ■ , u M - {v^}) = S^dua}; {v^}), (5.21) 

whereas the first one, 

G (o) 0'i. ■ • -M = 0i, ■ ■ ■ ,3m\B{v { ?) ■ ■ ■ B{v§)\ t) (5.22) 

is closely related to the partition function Z N (I5.16P of the Gaudin model. 

We then perform the summation in fl 5 . 1 9 [) and compute successively the functions G^ % \ 
Details are given in Appendix B. Finally, the scalar product S 1,1 ({u a }; {y\ }) has the fol- 
lowing determinant representation 

Sm(M, R W }) = G (M) K • • • ,u M ; {v^}) 

det^(K};fa (1) }) ( , „x 

n fc<i sin(w fc - Uj) sm(u k + w^) Y[ a>f3 sm{v K a - vy) sin(vi + v ( p ') 
where the matrix ^ /1 ' 1 ({m q ,}; {f} 1 }) is given by 

~i , m sin(2i> (1) ) sin(2ii a )F| 1) (M a ; {zA, {t> (1) }) 

sm(A 2 + £ - v) L) ) sin(Ai + £ - u} ; ) 
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Here we have introduced the functions 



(1)1 



sin(Ai + £ - u a ) sin(A 2 + £ - u a ) Ylf =1 sin(v^ - u a ) Y[™ =1 Bm.{v^ } + 



rM 



sin 2 (f - x) - u a ) sin 2 (u^ + u a ) 



1 - A 



+ 



1 + A 



sin(Ai + £ + u a ) sin(Ai + £ — sin(A 2 + £ + ««) sin(A 2 + £ - u a ) 

2M j M 

^ sin (7; . — 7-7 sin (7/ -I- 7, A ^ 



^ sin(7j a - z fc ) sin(u Q + z fc ) ^ sin ( Ma _ yW) sin ( Ma + v £h 



(5.25) 



Similarly, we obtain the scalar product S 2,2 ({u a }; {v^}) in terms of determinants with 



the parameters {v^ } satisfying the second set of Bethe ansatz equations f)4.10p . 

^ 2 ' 2 (K};{^ 2) }) 

det^KHtf}) 

Uk<j sin K - uj) sin(w fc + Uj) U a>( 3 sin(vi 2) - vf) sin(t»i 2) + vf) 

where the M x M matrix ^V 2,2 ({u a }', {^| 2 ' ) }) is given by 

_ sin(2tfj 2) ) sm(2u a )F^(u a ; { Zi }, {v^}) 



Here 



{*«) 



a,3 



sin(A 2 + £ + vf') sin(A a + £ + v) 



(2)- 



(5.26) 



(5.27) 



sin(Ai + ^ + n a ) sin(A 2 + £ + tt a ) Qjli s in K - M a ) llfcli sin K + M c 

sin 2 (i^ — u Q ) sin 2 (i^ + u a ) 
1 + A 1 - A 



+ 



sin(Ai + £ + u a ) sin(Ai + £ - u a ) sin(A 2 + £ + u a ) sin(A 2 + £ — u c 
2M m ^ 

+ ^sm(u a - z k ) sin(n Q + ^ sin ( Ma _ v ( 2 ) ) sin ( ?jQ + v ( 2 ) ; 



(5.2J 



6 Conclusions 

We have studied the XXZ Gaudin model with generic non-diagonal boundary terms. In 
addition to the inhomogeneous parameters {%}, the associated Gaudin operators {Hj}, 
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(13. 7p . depend on four free parameters {Ai, A2, £, A}. Thus our Gaudin operators are four- 
parameter ({Ai, A 2 , £, A}) generalizations of those in [29], three-parameter ({Ai, A 2 , £}) 
generalizations of those in [2U [3], and one-parameter ({A}) generalizations of those in [9]. 
The common eigenstates (Bethe states) of the operators are constructed by algebraic Bethe 
ansatz method. We have obtained the determinant representations (I5.16p . (I5.23|) and (I5.28j) 
of the scalar products for the boundary XXZ Gaudin model. 
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Appendix A: Derivation of (15.161) 



Take the partition function Z^\{v a }; {zi}) as an example, which satisfies the recursive 
relation, 

(2) ~ 
Z N ({M; {Zi}) = (Jl,"-,jN\B{v N )\ji,"-,jff-i) 

x ( 3l , ■ ■ -J^BivN^) ■ ■ • 5(^)1 it), (A.l) 

namely 

7 (2)/r X r X x v> sin(Ai + € + Zj) sin(A 2 + £ - z { ) sm(2v N ) 

^ sm(Ai + £ — v n) sm(A 2 + £ — v N ) sm{VN - z i) sin(t>jv + z%) 

x4 2 'i(KWw;{zjy. (A.2) 

The partition function Z N ({v a }; {^}), for any positive integer N, can be uniquely deter- 
mined by the initial condition: Z^ 2 \{v a }; {zi}) = 1 and the recursive relation (1A.2|) . 
Let us introduce a series of functions {Kj({v a }', {zi})\I — 1, • • • , N}, 

K N ({v a }\ {z^) 
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]la=i nil sin(w a + Zi) sin(« a - 



x det 

and prove the relation 

r(2) 



n„ >( g sin(w a - t^) sin(i; a + vp) Y[ k<j sin(z k - zj) sin(z k + Zj) 
sin(Ai + £ + 2j) sin(A 2 + £ - z,-) ■ 2u Q 



sin(Ai + £ - f a ) sin(A 2 + £ - v a ) sin (t> Q - Zj) sin (t> Q + z,,-) 



(A.3) 



z i ({ v a}] {zi}) = Kj({v a }; {z^), for any positive integer J. 



(A.4) 



We use the induction method. 



We have, for the case of N = 1, 



Z[ 2 \vi, zi)=.Ki(vi; z\)=— 



sin(Ai + £ + zi) sin(A 2 + £ — Zi) sin(2fi' 



sin(Ai+£— vi) sin(A 2 +£-vi) sia(vi—zi) sin^i+^i' 



(A.5) 



Suppose that f[Q|) holds for / < JV - 1. We prove that flA.3|) also holds for I = N 
holds. This can be done as follows. The determinant representation of K]s[({v a }; {zi}) 
implies that it satisfies the following recursive relation 

K N ({v a }; {zi}) 



N 

E 



x 



X 



sin(Ai + £ + z^ sin(A 2 + £ — sin(2t; A r) 
sin(Ai + £ - sin(A 2 + £ - Ujv) sin(t»jv - z { ) sm(v N + z { 

1 sinft^ - 2j) sin(f Z + Zj) -p-r sinfev - 0,-) shifty + zA 



n 



n 



Z;) sin(^ + ;& 



(A.6) 



The determinant representation (1A.3|) . the recursive relation (|A.6[) and the recursive 
relation (1A.2|) mean that -K0v({-i> a }; {^}) and Zj^({i> a }; as functions of vn, have 

the same set of simple poles, 



±z l ,X 1 + £,A 2 + £ mod(2vr), « = 1, • • - ,iV, 



at which both functions have the same residues. Moreover we can show that 



(A.7) 



4 2) (K};{*})I 



= ^(K};{zt})| 



(A.8) 



We thus conclude that ()A.4j) also holds for I = N. This completes the induction. 
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Finally we get the determinant representation of the partition function Z N {{v a }\ {zi}), 

7 Q)(S„\.S*X\- Illi Ilili sinfoq + Zi) sin(v a - Zi )det,yK^ 

[[ a>j3 sm(v a - vp)sm(v a + vp)[[ k<j sm(z k - zj) sm(z k + Zj) 

where the N x N matrix ^^ 2 \{v a }; {zj}) is given by 

^HM; = . ^^;;» + {p)Bm(2^) (A. 10) 

sm(Ai+^— f a ) sm(A 2 +f— ■u a ) sin (v a -%)sm [v a +Zj) 



Appendix B: Proof of (15.231) 



We prove (I5.23j) by calculating the functions ( I5.19P recursively. 

We illustrate our derivations for We firstly express in terms of G<-°\ 

G^( Ul \j 2 ,---,j M ;{v^}) 



0a, • • • >iM|C'(«i)|i, J2, ■■■,3m) x G {0) (j,j 2 , ■ ■ -J M ; {v^}) 

sin(Ai + £ — ^) sin(A 2 + £ + sin(2ui) 
. , ^ . sin(Ai + £ + sin(A 2 + £ + «i) sinfwi - sinfui + ^) 

x Uti ngi sin^ + g) sin^ - gQdet^j^ W) (R1) 
U a >p sintf } - uj ) sin(wi 1} + n fc<j sin(2 fc - Zj) sin(z k + Zj ) ' 

where ^({v^}; { Zi }) is a M x M matrix with matrix elements ^V({va ^}; {^i})a,j given by 



^');M)^ ,„ f ( ^^r' ) f ( ^:w 3in( ^l,., <M 

sm(Ai+£— Va ') sm(A 2 +£— fa ) sin [Va —Zij&m (f« ' + Zi 

Further, performing the summation and absorbing the results into the element of the first 
column of the matrix JV^ 1 ' below, we have 

G^imlh,-..,^^}) 

Ug=i Uk=2 sin(^ 1} + z ik ) sin(^j 1} - ^Jdet^WQ^}; u u i 2 , • • • , i m ) ^ 
llAf> a >/3>i sin(z;i 1) -^ 1) ) sin(wi 1) +^ 1) ) n 2 <jfc<i<M sin (^ -«0 sin (^+^) ' 

where 

^ = -7 m ; ^ for b > 2, (B.4) 

sm(ui - ZiJ sm(Mi + 
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M 



(1) 



Y]_ sin(wi - Zi k ) sin(ui + z ife 



fc=2 



M 



X 



X 



X 



z h ) sm(v^ - z h ] 



n 



M . / _ 



sin(z ifc + z 



sin(Ai + £ + £jj sin(A 2 + £ - 2^) sin(2^ 
sin(Ai+£--ui 1) ) sin(A 2 +£-wi 1) ) sin^fi^-ZjJ sin^ui^+ZjJ 
sin(Ai + £ - ZjJ sin(A 2 + £ + 2jJ sin(2ui) 



for 6 = 1. (B.5) 



sin(Ai+£+Mi) sin(A 2 +£+Mi) sin(wi— z h ) sin^i+^J 

Some remarks are in order. The summation in (1B.5|) is not over the full set of values, 
namely, i\ =fi i 2 ■ ■ ■ iu, (cf. [30J ) . Also there are second order poles (see ( 1B.6|) below), which 
make the calculation tedious. 

Taking advantage of the analytical properties of the trigonometric function and keeping 
in mind that {v!- } are solutions of the first set of Bethe ansatz equations, we can calculate 
the sum in (El to get 



M 

E 



n^LiSin(^ lj + gjjsin 



(1) 



n fc >i sin^! - Z(J sin(2; ifc + z iiy 

sin(Ai +£ — z il ) sin(A 2 + £ + z ix ) sin(2u 1 ) 
sin(Ai + i + ui) sin(A 2 + £ + «i) sin(«i - ZjJ sin(ui + ZjJ 

sin(Ai + £ + ZjJ sin(A 2 + £ - z^) sin(2va ) 
sin(Ai + £ - f o ) sin(A 2 + £ - t4 ) sin 2 (t>i - z^) sin 2 (t;i 1) + z^' 
^Kijui, {va ) }; {z h }) 

Uk=2 sill ( M l - Z ik) sin ( M l + ZiJ 



+ J2o tb ^£\u 1 ,{v^};{z h }), (B.6) 



b=2 



with 



sin(Ai+£-Mi)sin(A 2 +£-Mi)n^iSin(vi 1) -Mi)n*liSin(vi 1) +ui)sin(2vi 1) )sin(2ui) 



sin(Ai + £ - ) sin(A 2 + £ - Va 1,1 ) sm 2 ^ 11 - «i) sin 2 (ua AJ + Wi) 

2Af „ M 



1 



E- 



,2/.,(l) 
'a 

2 



y 

^ sin( Ul - z k ) sin(«i + 2 fc ) ^ s i n ( Ul _ y«) sin ( Ul + ^ 
1 - A 1 + A 



(B.7) 



sin(Ai + £ + Mi) sin(Ai + £ — Mi) sin(A 2 + £ + wi) sin(A 2 + £ — u\) _ 

where a& are coefficients which do not depend on a. Then only the first term of the r.h.s. of 
( 1B.6|) has nonzero contribution to the determinant of jV . 
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To show ( 1B.6j) . we first prove the following relation: 



M 

E 



n^Li sin(t;^ J + z h ) sin(t> 

r M 



(1) - z 



i^i 2 ,..,i M U k >i sin (^i - z i k ) sin (^ fc + z 

sin(Ai + £ - z h ) sin(A 2 + £ + ZjJ sin(2wi) 



x 



X 



sin(Ai + £ + ui) sin(A 2 + £ + «i) sin(-ui - z^) sm(ux + ZjJ 

sin(Ai + £ + ZjJ sin(A 2 + £ - gn) sin(2t;i 1) ) 

sin(Ai + £ - Ua ) sin(A 2 + £ - Ua ) sin 2 (wi 1) - ZjJ sin 2 (fi 1) + ^ 
^isin(tti, {v™}; {^}) 1 



nfl 2 Bin(«i - z lk ) sin(ui + sin(Ai + £ + ui) sin(A 2 + £ + m) 



y— — 

f * Gin 

+E 



.6=2 



sin (m - z ib ) Y\ k=2 ^ b sm{z ib - z ik ) U k =2 sin ( z h + z h) 

i I ^!i ) («i,tf > };W) 



6=2 



sin(«i - z ib ) I sin(Ai + £ + z ik ) sin(A 2 + £ - z. 



9 /riaii sin ( v Q 1)_M i) sin(vi 1) +-u 1 ) sin(2«i)nL 1 sin(A J -+£+«i) sin(A j +£-u 1 ) > 



9«i 



sin(2z i6 ) 



M 

' » Bin I 



u 1 =z i 



6=2 

' 2M 



sin( Wl - z*J H A k I =2 k ^ b sm(z lb - z ik ) ]jf =2 sm(z ib + z ik ] 
I " 2 



x 



^— ' <si 



+ 



^ sin(^ b - z k ) s\n(z ib + z k ) 
l - A 



E 



fe =i sm(z it -^')siii(2 i6 + v k 



+ 



l + A 



sin(Ai + £ + ZjJ sin(Ai + £ - sin(A 2 + £ + sin(A 2 + £ - z ib ] 

+ trigonometric polynomials with poles of sin(iti + z ik ), 



(B.8) 



where R(zi b ) is 



i2(; 



M M 

= sin(Ai + £ - z i6 ) sin(A 2 + £ + ZjJ JJ sin(v^ l) - ZjJ JJ sin(^ l) + z ib ) sin(2zjj 
sin(Ai + £ + z ib ) sin(A 2 + £ - z ib ) sin(2fj 1) ) 



?in(Ai + £ + z ib ) sin(A 2 + £ - z ib ) sin(2t;i 1) ) 

(Ai + £ - sin(A 2 + £ - u^) sin 2 ^ - z ib ) sin 2 (v£ ] + z, 

M M 



= sin(Ai + £ - ZiJ sin(A 2 + £ + ZjJ sin(v (1 - i 



Af 

Q=l 

(B.9) 



Equation( 1B.8l) can be proven by considering both sides as meromorphic functions of u±. 
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Then both sides of (1B.8j) have the same set of simple poles 



1 



ui = ±z il , mod(27r), where i\ ^ %i, • • • , Im, 

at which both have the same residues. Moreover, both sides tend as U\ — > oo. We thus 
conclude that relation (1B.8j) holds. Note that that except the first one on the RHS of (IB. 81) 
all other terms can be expressed in the form of Ylb a i>° / Kfb \ U U { v <* Y> { z i})- We thus obtain 

(El]). 

Repeating the above steps, by means of the recursive relations and changing the form of 
the determinant column by column, we obtain the final expression for S({u a }', {^>i })• 
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